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Abstract. We apply the Levenberg-Marquardt minimization algorithm to seismic and classical observables of the o-Cen binary 
system in order to derive the fundamental parameters of orCenA+B, and to analyze the dependence of these parameters on 
the chosen o bservables, on their uncertainty, and o n the physics u sed in stellar modelling. The seismological data are those by 
iBouchv & Carried feOO 3) for aCen A, and those bv lCarrier &Bourbaijl2003l) for aCen B. We show that while the fundamental 
stellar parameters do not depend on the treatment of convection adopted (Mixing Length Theory - MLT - or "Full Spectrum of 
Turbulence" - FST), the age of the system depends on the inclusion of gravitational settling, and is deeply biased by the small 
frequency separation of component B. We try to answer the question of the universality of the mixing length parameter, and we 
find a statistically reliable dependence of the or-parameter on the HR diagram loca tion (with a trend sim ilar to the one predicted 
by iLudwig et all 19991) . We propose the frequency separation ratios bv lRoxburgh & Vorontsovl 120031) as better observables to 
determine the fundamental stellar parameters, and to use the large frequency separation and frequencies to extract information 
about the stellar structure. The effects of diffusion, and equation of state on the oscillation frequencies are also studied, but 
present seismic data do not allow their detection. 
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1. Introduction 

a Cen AB is the binary system closest to the Earth (d=1.34pc). 
It sho ws an eccentric orbit ( e-0.5 19) with a period of almost 80 
years ( Pourb aix et al 1200*2 ). aCen A is a G2V star and aCen B 
a K1V one, sligthly hotter and cooler respectively than the 
Sun. Thanks to the high apparent brightness (Va = -0.01 and 
Vb = 1.33) and to the large parallax of the components, their 
stellar parameters are among the best known of any star ex- 
cept the Sun. The binarity, their well determined characteristics 
and the solar-like oscillations detected in both stars, provide a 
unique opportunity to test our knowledge on stellar evolution 
in conditions slightly different from the solar one. 

As a consequence, a great number of theoretical stud- 
ies dealing with aC e n has been published since the one 
by iFlannerv & Avresl Jl978l) (see lEggenberger et al] |2004) 
for a comprehensive review). Before the definitive identifi- 
cation of p-mode frequenci e s in t he aCen A power spec- 
trum by Bou chv & Ca rrier (2002J), the uncertainty in the 
parallax (and therefore in the masses) and in the chemi- 
cal composition did not allow an unambiguous determina- 
tion of stellar parameters. Some controversial results came 
up concerning, for instance, the universality or not of the 
mixing-length parameter (or) describing the stellar convection 
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dNoels et alJ Il99l| Edmonds et alJ fl99l Lvdon et al.1 Il993l 
i Neuforgd 1 19931 Morel et al l 12 000: TFCTnandeT^TNeuforga 
1199.4 iGuenther A Demanm d bOOd) : the role of the chemi- 
cal composition on discriminating between these two possi- 
bilities IFernandes & Neuforge 1995), and on the presence or 
not o f a convective core, and its e ffect on the ag e of the sys- 
tem (Guenther & Demarque 2000). S ome efforts jBrown et alJ 
Il994l lGuenther & Demarquel l200ri iMorel et alJ l2000h were 
also devoted to st udy the capability of solar- like oscillations 
expected in aCen ( Kieldsen & Beddi ng! 1995b to constrain the 
fundamental stellar parameters and the physics included in stel- 
lar models. 

In addition to th e p-mode identification by 
IBouchv & Carried d2002l) . IPourbaix et ail J2002I) improved 
the precision of t he orbital paramet ers and, adopting the 
parallax derived by Soderhielm ( 1999), provided very precise 
masses for orCen A and B. These high quality data stimulated 
new calibra t ions o f the system by Theveni n et alJ (2002) and 
iThoulet al] (|200 3|). The two teams reached different results. 
While lThe\^nrn et al ] ll2002l) could not fit the seismic data 
without changing the masses more than 4cr with respect to 
Pourbaix's data, the second group fitted the cCen A p-mode 
spectrum and the spectroscopic constraints using a single value 
of mixing length parameter, and with the new values for the 
masses. 
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Interferometric measurements with VINCI/VLTI by 
iKervella et all J2003I) have provided high precision values of 
the angular diameter of aCen A and B, and new observations 
bv lCarrier & Bourbanl (l2003) have allowed to identify p-mode 
frequencies also i n the B component. The se new constraints 
have been used bvfeggenber ger et"al] J2004I) . Their calibration, 
based on a grid of models obtained by varying the mixing- 
length parameters, the chemical composition, and the age, 
leads to a stellar model in good agreement with the astrometric, 
photometric, spectroscopic and asteroseismic data, and they 
assert that "the global parameters of the aCen system are now 
firmly constrained to an age of t — 6.52 + 0.30 Gyr, an initial 
helium mass fraction Yj — 0.275 + 0.010 and an initial metal- 
licity (Z/X)i = 0.0434 ± 0.0020" and that "the mixing length 
parameter a of the B component is larger than the one of 
the A component" . Th ese results are quan titatively consistent 
with those obtained bv lThoul et al.l (Eo03). nevertheless, both 
groups have performed the calibration assuming a given, but 
different, "physics". 

The aim of this work is to study, following the theoretical 
analysis of the utili ty of seismology to constrain fundamental 
parameters made bv lBrown et al.ldl994l) . the dependence of the 
set of parameters obtained by fitting the observables, on the 
details of the fitting procedure. That is: i) the kind of constraints 
included in the x 1 functional that drives the fitting procedure, 
as well as other effects of their uncertainties, and ii) the physics 
included in the stellar evolution theory. Only in this way we 
will be able to provide an estimation of the uncertainty in the 
obtained set of stellar parameters, and to evaluate the degree at 
which the present data precision can constrain the physics in 
stellar evolution models. 

To do that, and also with the prospect of dealing with a great 
quan tity of seis mic data from sp atial missions such as M OST 
jMatthewslll998i) and COROT (Bagl in & The COROT Teaml 
199£ jj) and from new generation spectrograph such as HARPS 
(Bouchv 2002), we have implemented a non-linear fitting algo- 
rithm that performs a simultaneous least-square adjustment of 
all the observable characteristics, the classical and the seismic 
features. 

The fitting method is described in Sect. |2] In Sect. [3] we 
discuss the different sets of classical and seismic observables 
that will be included in our^ 2 quality function. The physics in- 
cluded in the stellar evolutionary code is summarized in Sect.0] 
The results of these different combinations of observables and 
of different physics are presented and discussed in Sect. |5] A 
special effort has been devoted to the problem of stellar con- 
vection (Sect l5.2t . We will try to answer to the question about 
the universality of the mixing-length parameter, and to study 
the effect of different convection treatments. With respect to 
the convection modell ing, a controversial result was obtained 
bv [Morel et al . (2000), they reached different ages depending 
on whether they used the classical MLT theory ( Bohm- Vitense 
Il958l) . or the FST theory bv lCanuto fc MazzitelliMl99lUl992l 
hereafter CM91,CM92). Hence, we performed different cali- 
brations changing the convection treatment (FST or MLT), as 
well as different MLT calibrations either with a unique or dif- 
ferent a values for each component. The effects of using dif- 
ferent equation of state, including or not gravitational settling, 



and adopting a different solar mixture are discussed in Sect. 15.31 
Finally, results and conclusions are summarized in Sect. [6] 



2. Calibration method 

Usually the approach to analyze stellar oscillation data is rather 
conventional: i) several stellar models that bracket the known 
observational constraints (typically composition, mass, lumi- 
nosity, and effective temperature) are computed ; ii) p-mode os- 
cillation frequencies are calculated for the models; and Hi) the 
model oscillation spectra are compared with the observed one. 
We believe that a different approach is needed, so that astero- 
seismology is directly included in the calibration procedure and 
the results are not biased by a limited/subjective exploration of 
the parameter space or strongly depent on an initial guess of 
the model parameters. 

The development and use of objective and efficient proce- 
dures to fit stellar models to observations has become of evident 
utility in particular since seismic constr aints are included in the 
modelling (see e.g. iBrown et alJl994l) . 

Guenther & Br ownl J2004I) have proposed a method that 
quantifies at which degree the oscillation spectra as obtained 
from a grid of models parametrised in mass, age and composi- 
tion reproduce the observations. Models providing a minimum 
in the^- 2 , defined by the differences between the theoretical and 
observational frequencies, are s elected for an additional inspec- 
tion in a finer grid. As noted bv lGuenther & Brownl d2004l) the 
first problem in this kind of approach is its computational cost. 
Moreover the direct fit of the oscillation spectra implies a good 
knowledge of the surface layers o f the star, as that s t rongly 
affects the exact frequency values. Guent her & Brownl (12004) 
quantify the uncertainty in the theoretical frequencies due to 
our poor modelling of the external layers by using the dis- 
crepancies between observed and theoretical solar frequencies. 
But, how good is this estimation for stars with different superfi- 
cial gravity, chemical co mposition and age? On the other hand, 
Eggenberger et alJ d2004l) use the aforementioned conventional 
approach, and only in a second and a third phase take into ac- 
count the asteroseismic data, first the large and small frequency 
separations (Av, Sv), and then the frequencies. 

Here we propose a calibration method that finds the pa- 
rameters of the system by the minimization of a y 1 functional 
including at the same time classical and asteroseismic observ- 
ables. The parameters of the models are, as usual, the mass, 
initial chemical composition, age, parameter(s) of convection 
a. The observables could be chosen among the masses, T e ff, L, 
R, [Fe/H], Av, 6v (or combinations of these frequency differ- 
ences). Of course, being o-Cen a binary system, the same initial 
chemical composition and age have to be assumed when cali- 
brating components A and B. 

We define a quality function measuring the distance be- 
tween models and observations, that is, a goodness-of-fit mea- 
surement by: 
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where O , cr and O are respectively the observed value, 
observational uncertainty and theoretical prediction of each of 
the N D (A+B components) observables considered. 

The choice of the observables included in the objective 
function is thoroughly described in Section[5] 

In the general case of a binary system the model that gen- 
erates the observables and their derivatives has seven free pa- 
rameters (or six if qta = cb)- The most substantial part of the 
model consists of a stellar evolution code (CLES, see Sec. 
which takes as inputs the masses of each component (Ma, Mb), 
the initial chemical composition of the system (Y, Z), the age 
and the convection mixing-length parameters {a a, cb), that in 
general are assumed to be different for the two stars. 

We compute oscillation frequencies for the models by solv- 
ing the equations of adiabatic oscillations (OSC), and deter- 
mine (Av) and (6v) for the degrees t — 0, 1,2, 3. This is not 
done by a least square fit to the computed frequencies, based on 
the asymptotic properties of low degree modes, but by making 
the average of the theoretical separations in the domain of ob- 
served radial order n(n = 15-25, for aCen A, and n = 17-27 
for aCen B). The observed average separations have been de- 
termined from observed frequencies in the same way. 

For each component the evolution code provides the stel- 
lar luminosity, the radius and the effective temperature, as well 
as the model quantities required in the subsequent calculations 
of the oscillation frequencies. In the calibration process, the 
derivatives of the observable quantities are obtained varying 
each of the parameters (Ma, Mb, Z, Y, a a, ^b, and the age). 
We do not derive colors or visual magnitudes, and we do not 
include orbital elements such as apparent semi-major axis (a'), 
orbital period (P nr h) or parallax . We assume the parallax as de- 
termined by Soderhielm ( 1999) and the values of observables 
based on it. The masses, however, are considered as parameters 
and also as observable quantities in our calibrations. 

In the calibration of a binary system the large number of 
variables involved, both in terms of model parameters and ob- 
servables, suggests the use of a least-squares based fitting pro- 
cedure. This is particularly useful in this kind of calibration as 
we fit at the same time classical and seismological observables 
without making first a selection based on the HR location of 
the system. 

As shown in lBrown et alJ i ll 99 4) the observable quantities 
depend on the parameters in a complex way: stellar evolution 
is not a linear problem. Most of the observables are influenced 
by several parameters, and hence the connection between ob- 
servables and parameters that could conceptually seem to us 
the simplest one will not always provide the correct results. 

2.1. Optimization algorithm 

We use the gradient-expansion algorithm known as Levenberg- 
Marquardt method. This algorithm combines the advantages of 
an expansion method, i.e. rapid convergence close to the min- 
ima, with those of the gradient-search, that is, a rapid approach 
to a far away minimum. This method has as well the strong ad- 
vantage of being reasonably insensitive to the starting values of 
the parameters. 



At each step the fitting function is linearized calculating 
numerical derivatives (centered differences) of the observables 
with respect to each model parameter. The displacement in the 
parameter space, leading to a lower value of y 2 , is calcu lated 
following the prescription of Bevington & Robinson (2003), p. 
161-164 and the iterative procedure is ended when^ 2 no longer 
changes more than 2%. In calibrations with Ma,Mb, Z, Y, cta, 
as, and the age as model parameters, convergence is typically 
achieved in 3-4 iterations and at each of them the computation 
of 16 evolutionary tracks is needed to evaluate centered deriva- 
tives. The result of such a local minimization could be sensitive 
to the initial guess of the parameters; therefore, in order to get 
a more reliable final solution we perform several runs starting 
from different points in the parameter space. The effect we find 
is limited to a variation of the number of iterations needed to 
reach the minimum ^ , whereas the final parameters of the sys- 
tem differ much less than their uncertainty. On the other hand, 
building a 6-dimensional dense grid of models seems impracti- 
cal considering the aim of evaluating the effects of using differ- 
ent physical prescriptions in our models (e.g. different equation 
of state, different metal mixture etc). 

It is sometimes nonetheless possible to make fairly direct 
connections between the observables and the parameters, par- 
ticularly when one observable is much better determined than 
the rest. The solution is determined by the parameter that is 
known with very small uncertainty. The uncertainties in the pa- 
rameters for these fits are calculated from the diagonal terms in 
the error matrix (inverse of curvature matrix in the parameter 
space) and are, in general considerably larger than the uncer- 
tainties obtained in the grid- and gradient-search methods. The 
latter are obtained by finding the change in each parameter to 
produce as change of x 2 of 1 from the minimum values, with- 
out re-optimizing the fit, while there is a strong suggestion that 
correlations a mong the parameters play an important role in 
fitting (see e.g. Bevington & Robinson 2003). 

3. The choice of the observational constraints 

3.1. Non-asteroseismic constraints 

Due to the proximity of the cCen binary system, the preci- 
sion on the measurement of its trigonometric parallax is po- 
tentially very high. Unfortunately, some discrepancies have ap- 
p eared among the mos t recent published values (see Table 10 
in lKervella etaIll2003b . lGuenther & Demaraud ( 120001) studied 
the uncertainty in the stellar parameters due to the different par- 
allax values. This, indeed affects the determination o f mass , lu- 
minosity and radius. Following Eggenb erger et alJ (2004) we 
adopt 7i = 747.1 ±1.2 mas (Soderhielm 1999) and, therefore, 
the c orresponding mass values determined by IPourbaix et alJ 
(120021) : (M A = 1.105 + 0.007 M Q , M B = 0.934 + 0.006 M Q ) 
and the radii : (Ra = 1 .224 + 0.003 R Q ; R B = 0.863 + 0.005 R Q ) 
dKervellaetal.l2003h. 

As in the case of the parallax, there is a large scatter in the 
published values of other quantities, such as T e ff, luminosity, 
and metallicit y. We d e cided to use the same values adopted 
by Eggenbergeret al. (2004) in order to have a reference 
model. Eggenberg er et al.l IJ2004) took as T e g for the compo- 
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nent A a value to encompass t hose given by two spec troscopic 
deter minations, the one from INeuforge-Verheecke & Magain 
J 19971) (T.« ^=583 0+30 K, r pffR =5255+50 K), and the one by 
Mor el et alT (l2000) based on a re-analysis o flChmielewski et"all 
i 19921) spectra (r eff A =5790±30 K, r pff R =5260+50 K), use d re- 
spectively in the aCen calibrations by Th oul et alJ {2003) and 
IThevenin et a J2002). We will use the effective temperature 
as constraint in our minimization method only in one of the 
calibrations (Alt,Blt, in Table |2j, since the precise determi- 
nation of the radius provides a narrower domain in the space 
of observable quantities . The luminosity values adopted by 
Egg enberger et alJ J2004I) come from a new and weighted cali- 
bration of previous Geneva photometric data, where they have 
also coherently taken into account the effective temperatures 
an d the parallax. These v alues cover the domain considered 
by IThevenin et alJ J2002I) that considered an error bar twice 
smaller, and a lower luminosity for the component B. The lu- 
minosity values, directly determined from the adopted radius 
and effective temperatur e, are in very good agre ement with 
the those determined by Eggenber ger et alJ (2004): L A /L G = 
1.518 ± 0.06, Lb/L = 0.507 + 0.0 25. On the other hand, the 
values determined by Piipersld2003h for La are much larger and 
only marginally overlap the values considered here. 

The precise values of masses and radius provide also pre- 
cise values of the surface gravity for both stars: log gA=4.305 ± 
0.005 and log g R = 4.536 + 0.008, while the spectroscopic val- 
ues determined by iNeufor ge-Verheecke & Magainl Jl997h are 
loggA = 4.3 4 + 0.05 and logg R = 4.51 + 0.08. These values 
were used by Tho ul et alJ (12003) to fix the luminosity domain 
, leading to higher central valu es and to larger error bars with 
respect to those deter mined bv [Eggenber ger et al.1 d2004l) and 
IThevenin etail d2002l) . 

Also for the metallicity of both components there is no 
complete agreement in the literature: [Fe/ HU = 0.20 + 0.02, 
and [Fe/H] B = 0.23 ± 0.03 from iMorel et all feOOOl) . and 
rFe/Hl a = 0.25 ± 0.02, and rFe/Hl R = 0.24 + 0.03 from 
INeuforge-Verheecke & Magainl dl997l) . The uncertainty in the 
observable Z/X is quite large, if we take into account also 
the 10% in the (Z /X)^. We have taken the value adopted by 
iThoul et alJ l200A that is Z/X = 0.039 ± 0.006, the same for 
both stars. T he detailed abundan ce analysis of a Ce n A and B 
earned out in lNeuforge-Verheecke & Magainl d 19971) suggested 
no evidence for a different metal mixture relative to the sun, 
therefore a ll our models were comp uted assuming the solar 
mixture bv iGrevesse & Noels! (^993), except for the calibra- 
tion (A5,B5) in which we have considered the recently d eter- 
mined solar metal abundances ( Asplund et al. 2004, 2005) that 
implies (Z/X) Q = 0.0177. 

The non-asteroseismic constraints used in this work are 
summarized in Tabled 

3.2. Seismic constraints 

Solar-like oscillations generate periodic motions of the stellar 
surface with periods in the range of 3-30 min and with ex- 
tremely small amplitudes. Frequency and amplitude of each 
oscillation mode depend on the physical condition prevailing 



Table 1. Non-aster oseis mic constr aints. 

Referencesm lEggenberger et al] 120041) : f2i: lThoul etaiN2003h . 





A 


B 


Ref 


M/M 


1.105±0.007 


0.934+0.006 


(1) 


T eg 


5810+50 


5260+50 


(1) 


R/R 


1.224 ± 0.003 


0.863+0.005 


(1) 


L/L 


1.522+0.030 


0.503+0.020 


(1) 


Z/X 


0.039+0.006 


0.039+0.006 


(2) 



in the layers crossed by the waves and provide a powerful seis- 
mological tool. Helioseismology led to major improvements in 
the knowledge of Sun structure and to revision of the "stan- 
dard solar model". The potential utility of seismology applied 
to other stars, in particular arCen, to constrain the stel lar param- 
eters was extensively studied in Brown et ID C994), and also 
bv lGuenther & Demarqud ll2000). 

Several groups had made thorough attempts to detect the 
signature of p-mode oscillations in aCen A, but their result s 
were not confirmed. Only recently B ouchv & Carried (J2002), 
from high precision radial velocity measurements with the 
CORALIE echelle spectrograph have yielded a clear detection 
of p-mode oscillation, and identified several modes between 
1800 and 2900 yuHz , and with an envelope amplitude of about 
31 cms -1 . Assuming that freque ncy mode s v„e satisfy the sim- 
plified asymptotic relation (Tassoul 1980): 

vt„«Av \n + - +e\-e(e+\)-^ (2) 

and assuming the parameter e near the solar one (1.5) they es- 
timated: Av = 105.5 + O.lyuHz, 5v Q2 = 5.6 ± 0.7 /kHz, e = 
1.40 ± 0.02 and identified 28 p-modes with degree £ = 0,1,2 
and order between n = 15 and n = 25. 

Notice that the given errors come from the autocorrelation 
algorithm, but we must keep in mind that their frequency res- 
olution is only 0.93 yuHz, and that they derive an uncertainty 
in the frequency determination equal to 0.46 /iHz. They also 
point out that an error of ±1 .3^/Hz could have been introduced 
at some identified mode frequency, that could explain the dis- 
persion of mode frequency around the asymptotic relation. In 
particular higher observational uncertainty could affect mainly 
the 1-2 modes that determine the value of 5vo2 for the lower 
an d higher frequency £yp2(n= 16 and 25). 

ICarrier & Bou rban (2003) have also detected solar-like os- 
cillations in the fainter component, aCen B. Only twelve fre- 
quencies, between 3000 and 4600 /iHz have been kept in 
the final list of identified p-modes, four of them with a de- 
tection level lower than 3<x, and it is recommended to take 
them with caution. As for component A the frequency res- 
olution is 0.93 yuHz. The large and small separations, deter- 
mined by autocorrelation of the asymptotic relation, are re- 
spectively Av = 161.1 ± 0.1 //Hz and 6v 02 = 8.7 + 0.8 yuHz. 
We must note here that the value derived for dv^i comes from 
only few p-modes. In fact, from their frequency table is only 
possible to obtain two va lues: 6vm.(n = 21) = lO. OjuHz and 
6vq2(i - 23) = 7.0(uHz. [Carrier & BourbarJ d2003) expect a 
rotational splitting ~ 0.3^Hz that, given the frequency resolu- 



A. Miglio and J. Montalban: Constraining fundamental stellar parameters using seismology 



5 



tion, could imply an increase of the uncertainty of frequencies 
for modes of degree £ — 1 and 1 = 2. 

Very recently, new o bservations of this system by 
iKieldsen & Beddind (12004 have confirmed the values deter- 
mined by B ouchv & Carried §002) concerning the frequency 
separations of component A. However, component B with this 
new more precise data shows Avo = 161.4/^Hz and (5vo2 = 
lO.ljt/Hz. Our computations have been done before these val- 
ues were available, therefore, we will not take them into ac- 
count in our calibration. Notice that, nevertheless, those values 
are anyway reached and imposed by the other observables used 
in some of our calibrations (A3,B3), see Sec. [5] 

How should these seismological observations be used to 
constrain our stellar models? The classical way is to use the 
large and small separations to characterize the power spectrum 
of solar-like oscill ations. The sta ndard asymptotic theory of 
stellar oscillations llTassou]i ri980) relates the averages values 
of high-radial order / low-degree small and large separations 
to conditions in the st ellar core (6v) and to the m ean densi ty of 
the sta r (Av). Recently|(juenther & Brown (2004) and Metcalfe 
(2005) have proposed to use directly the p-mode frequencies as 
observabl es to constrain the stellar models. 

Bro wn et"aD i 1994b theoretically analyzed the case in 
which the individual frequencies are included as observables. 
Their purpose was to illustrate the potential loss of information 
resulting from representing the spectra in terms of the large and 
small separations derived from the expected asymptotic behav- 
ior. The dominant source of frequency changes is very close to 
the stellar surface. It could be difficult to disentangle these ef- 
fects from the uncertainties in the treatment of the physics of 
the outer layers, where non-adiabaticity and dynamics effects 
of convection have to be taken into account. 

The oscillation frequencies, the large and small separations 
depend on the structure of both the inner and the outer lay- 
ers of a star, so model fitting and testing techniques to probe 
the interior structure of the stars are dependent on our having 
a good understanding of the structure of the outer layers. But 
these are just the layers where our ignorance is greatest; non- 
adiabatic convection is important but not understood, the oscil- 
lations are non-adiabatic in the surface layers and the structure 
of real stellar atmosphere is poorly understood. For example the 
oscill ation frequencies p redicted b y the re ference solar models 
(S96)(Christensen-Dalsgaar d et alJll99a) differ from the ob- 
served values up to IOjuHz at the higher end of the observed 
frequency range. 

In a first step, we will include as seismic constraints in our 
fitting algorithm the combinations of frequencies: the large 

Av„,f = V„,f - V n -\,t 

and small 

frequency separations defined from the identified p-mode fre- 
quencies for both co mponents. 

As discussed in | Christensen-Dalsgaard etafl Jl995t) and 
iDi Mauro et alJ J2003I) . care has to be taken when considering 
as a constraint in the modelling the large separation, as its av- 
eraged value at high frequencies could be influenced by near- 
surface effects as well. This is in fact the case when comparing 
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Fig. 1. Solar large frequency di fference Av„ i (upper panel) from 
standard seismic solar model S96 ( Christensen-Dalsaaard et al. 1996) 
(solid line), compared to the observational solar large separation (dots) 
iBasu et alll997l) . Lower panel: as upper panel but for the ratio ro2. 

the observed and the predicted low-degree large separations of 
the Sun (see Fig. 0, where the disagreement of the order of 
a //Hz is related to a simplified treatment of the model outer 
layers. 

With the aim of checking whether the calibration we per- 
form considering (Av) and (Sv) is not biased by a simplified 
treatment of the outer structure in our models, we considered 
the effect on the calibration of choosing as seismic constraint 
ro2, the ratio between the small and large frequency separations 
defined by: 

Sv„fi V„ n - Vr-1,2 



>"02 («) 



Av„ 



Vnl - Vn-1,1 



(3) 



for 6 different orders n of the compo nent A. This combinatio n 
of frequencies, as presented in Roxburgh & Vorontsov (2003), 
is to a great accuracy independent of the outer layers of the star, 
and therefore represents a reliable indicator of the conditions in 
the deep stellar interior. 

4. Stellar models 

All stellar model sequences are calculated using the CLES 
code (Code Liegeois d' Evolution Stellair e). The op acity ta- 
bles are those of OPAL96 (Igle sias & Rogers1ll99rj) comple- 
mented at T < 6000 K with lAlexander & Ferguson (1994) 
opacities. The relative mixture of heavier than helium ele- 
ments, used in the opacity and equation of state tables, is 
the solar one according to iGrevesse & Noelsl dl993h . The nu- 
clear e nergy g eneration rou t ines a re based on the cross sec- 
tions by Caughlan & Fowlerl dl988l) and screening factors from 
Salp etei Jl954l) . CL ES allows the choice between two equa- 
tions of state: CEFF ( Christ ensen-Dalsgaard & Dappenl ll992) 
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and OPAL01 jRogers & Navfonovl2002t) . Most of the models 
have been computed using OPAL01, but we have also made a 
calibration using CEFF to study the capability of seismology 
to constrain the EoS. 

All our stellar models were obtained from evolutionary 
tracks including the pre-main sequence phase, and ending at 
~9 Gyr. The stellar models have approximately 1200 shells, the 
last one corresponding to T=T e ff as determined u sing as bound- 
ary conditions those given by the iKuruczl ill 998h atmosphere 
models. Furthermore, for the computation of oscillations we 
have added atmospheric layers from T = T e ff up to r = 10~ 4 . 

We have computed models where the convection is treated 
both using the Mixing Length Theory (MLT iBohm-Vitensel 
Il958l) with the formalism described in ICox & GiuITi 19681) . 
and the Full Spectrum of Turbulence (FST Canuto et al. 1996) 
with a formal i sm si milar to the one used by Mo rel et alJ (J2000) 
or i Bernkopil Jl998h. T hat means convective fluxes as given 
by ICanuto et alJ ( 119961) . but another prescription of the scale 
length. The parameter a of the MLT and the corresponding in 
the FST are considered as free parameters of the model cali- 
bration, and the values obtained are compared with the values 
required in the solar calibration for the same physics. 

Finally we have computed stellar models with and without 
gravitational settling of helium and metal s. The micro scopic 
diffusi on formulation is that given bv lThoul et alJ i 19941) . solv- 
ing the Burgers ( 1969) equations for H, He and Z and thus con- 
sidering diffusion due both to thermal and concentration gradi- 
ents. We also assume complete ionization and the effect of ra- 
diative acceleration is ignored. This assumption is completely 
justified for the precision of t he models and for the masses con- 
sidered (Tur£2MitiiIi^98). 



0.2 



0.19 



5? 0.18 



0.17 



-0.26 



-0.28 



-0.3 - 



-0.32 



A 



A2 
A3 



3.77 



A3f 
A3c 
A3a ( 
A4 



A 5 
A3e 

Any 



i.765 
log T„ 



3.76 



0.34 




log T e , 



5. Results 

In Table 13 we list all the fits of the binary system aCen. The 
calibrations differ in the choice of the observables included in 
the objective function (x 2 ) and in the used physics. The first 
column in Table [2] identifies the models, the second, third and 
fourth columns describe the physics; the fifth one resumes the 
characteristic of the seismic and non-seismic constraints used 
in the calibration process, while all the following columns give 
the values of the parameters (and the uncertainty in each of 
them) providing in each case the minimum^ 2 . The values of 
the observable quantities derived for each set of parameters is 
listed in Table [5] There, the last column;^ 2 is the the value of 
the "reduced"^ 2 , defined as: 



x 



(4) 



N a -N p 

were N t) is the number of observables and N p the number of 
model parameters. 

In the three following subsections we shall analyze the ef- 
fect of changing the non-seismic constraints, the seismic con- 
straints and finally the physics used in the models, that is, a 
different treatment of convection, different equation of state, 
different solar mixture, models including or not gravitational 
settling, and the effect of considering overshooting in our mod- 
els. 



Fig. 2. HR diagram location of models with stellar parameters from 
differing fits. Labels corresponds to those in Tables [5] and [3] Upper 
and lower panel correspond respectively to components A and B. The 
error boxes for r eff , log L/L Q , and radii correspond to lcr (solid line) 
and 2cr (dashed-line) 

Fig.|5]shows the HR location of each of these fitted models. 
We have indicated the error boxes in T e $ log L/L , and radius, 
corresponding to lcr (solid line), and 2<x (dashed-line). 

5.1. Effect of classic and seismic observable quantities 

A first calibration (Alr,Blr) was performed including in the 
X 1 function the luminosity, the radii and the actual (Z/X) s , 
as well as the average large and small frequency separations 
(computed as described in Sect. for each component. The 
model parameters are the initial chemical composition (Yq,Zq), 
the mixing length parameters (ota, a R ) and the age (t) . In thi s 
calibration we consider, following Eg genberger et alJ 



that the masses are perfectly determined, and we assume the 
mass of each c omponent to be fixed t o its observational central 
value, as given Pourba ix et alJ d2002l) . The parameters provid- 
ing the minimum^- 2 are: r = 6.8 +0.5 Gyr, o-a=2.00, a^-2.2A, 
Fo=0.276 and Zo=0.0325. We note that t hese results are in 
comp lete agreement with those obtained by Eg genberger et alJ 
(2004): r = 6.5 + 0.2 Gyr, Z = 0.0302, Y = 0.275, and 
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Table 2. Sets of parameters of fitted models. Meaning of numeric labels: 1 (fixed masses and Av, Sv as seismic constraints); 2 (as 1 but with 
variable masses); 3 (as 2 but using ro2(n) as seismic constraint); 4 (as 3 but including convective overshooting) and 5 (as 3 but using Asplund et 
al. (2005) instead of Grevese & Noels (1993)). Meaning of alphabetic labels: nd (non diffusion models); f (FST convection treatment); e (CEFF 
EoS instead of OPAL01 one); c (a unique mixing-length parameter); ns (fit without seismic constraints); only for the case 1, t refers to effective 
temperature as constraint, and r to radius as constraint. 



Model 


Conv 


Diff 


EoS 


fitting 


M/M 


e 




e 


Age 


6 


Y 


e 


Zo 




Alt 


MLT 


Y 


OPAL 


r eff , M flx 


1.105 




1.99 


0.11 


7.0 


0.5 


0.273 


0.016 


0.0316 


0.003 


Bit 


MLT 


Y 


OPAL 


(Av), <<5v> 


0.934 




2.22 


0.07 


7.0 


0.5 


0.273 


0.016 


0.0316 


0.003 


Air 


MLT 


Y 


OPAL 


M flx 


1.105 




2.00 


0.11 


6.8 


0.5 


0.276 


0.015 


0.0325 


0.003 


Blr 


MLT 


Y 


OPAL 


(Av), (Sv) 


0.934 




2.24 


0.07 


6.8 


0.5 


0.276 


0.015 


0.0325 


0.003 


Alnd 


MLT 


N 


OPAL 


M 6x 


1.105 




1.84 


0.10 


7.1 


0.6 


0.265 


0.015 


0.0284 


0.003 


Bind 


MLT 


N 


OPAL 


(Av), (Sv) 


0.934 




1.99 


0.05 


7.1 


0.6 


0.265 


0.015 


0.0284 


0.003 


Ale 


MLT 


Y 


CEFF 


M flx 


1.105 




1.93 


0.10 


6.7 


0.6 


0.281 


0.015 


0.0324 


0.003 


Ble 


MLT 


Y 


CEFF 


(Av), (Sv) 


0.934 




2.08 


0.06 


6.7 


0.6 


0.281 


0.015 


0.0324 


0.003 


Al-3 


MLT 


Y 


OPAL 


M flx , 


1.105 




1.86 


0.06 


6.0 


0.4 


0.277 


0.017 


0.0302 


0.003 


Bl-3 


MLT 


Y 


OPAL 


r(n(n) 


0.934 




2.18 


0.05 


6.0 


0.4 


0.277 


0.017 


0.0302 


0.003 


A2 


MLT 


Y 


OPAL 


M var 


1.104 


0.006 


1.96 


0.11 


6.4 


0.6 


0.282 


0.016 


0.0326 


0.003 


B2 


MLT 


Y 


OPAL 


(Av), <<5v> 


0.926 


0.005 


2.11 


0.10 


6.4 


0.6 


0.282 


0.016 


0.0326 


0.003 


A2c 


MLT 


Y 


OPAL 


M val a A = a B 


1.107 


0.005 


1.78 


0.08 


5.4 


0.5 


0.284 


0.016 


0.0305 


0.003 


B2c 


MLT 


Y 


OPAL 


(Av), <<5v> 


0.919 


0.004 






5.4 


0.5 


0.284 


0.016 


0.0305 


0.003 


A2f 


FST 


Y 


OPAL 


Af var 


1.105 


0.006 


0.77 


0.05 


6.4 


0.6 


0.282 


0.016 


0.0325 


0.003 


B2f 


FST 


Y 


OPAL 


(Av), <<5v> 


0.927 


0.005 


0.89 


0.05 


6.4 


0.6 


0.282 


0.016 


0.0325 


0.003 


A3 


MLT 


Y 


OPAL 


M var 


1.110 


0.006 


1.90 


0.07 


5.8 


0.2 


0.284 


0.014 


0.0322 


0.003 


B3 


MLT 


Y 


OPAL 


ro2.(n) 


0.927 


0.005 


2.05 


0.08 


5.8 


0.2 


0.284 


0.014 


0.0322 


0.003 


A3c 


MLT 


Y 


OPAL 


M mr , a A = a B 


1.113 


0.006 


1.95 


0.05 


5.8 


0.5 


0.276 


0.017 


0.0317 


0.003 


B3c 


MLT 


Y 


OPAL 


min) 


0.922 


0.004 






5.8 


0.5 


0.276 


0.017 


0.0317 


0.003 


A3ff G 


MLT 


Y 


OPAL 


M var , a = a Q 


1.114 


0.006 


1.91 




5.7 


0.2 


0.283 


0.010 


0.0314 


0.002 


B3ff 


MLT 


Y 


OPAL 




0.921 


0.004 






5.7 


0.2 


0.283 


0.010 


0.0314 


0.002 


A3f 


FST 


Y 


OPAL 


M var 


1.110 


0.006 


0.74 


0.03 


5.8 


0.2 


0.285 


0.015 


0.0323 


0.003 


B3f 


FST 


Y 


OPAL 


rca(n) 


0.927 


0.005 


0.86 


0.03 


5.8 


0.2 


0.285 


0.015 


0.0323 


0.003 


A3nd 


MLT 


N 


OPAL 


M var 


1.108 


0.006 


1.84 


0.07 


6.3 


0.4 


0.270 


0.015 


0.0281 


0.003 


B3nd 


MLT 


N 


OPAL 


f02(n) 


0.929 


0.005 


1.99 


0.07 


6.3 


0.4 


0.270 


0.015 


0.0281 


0.003 


A3e 


MLT 


Y 


CEFF 


M var 


1.109 


0.006 


1.85 


0.07 


5.7 


0.3 


0.288 


0.015 


0.0322 


0.003 


B3e 


MLT 


Y 


CEFF 




0.927 


0.005 


1.93 


0.07 


5.7 


0.3 


0.288 


0.015 


0.0322 


0.003 


A4 


OV 


Y 


OPAL 


M var 


1.112 


0.006 


1.77 


0.05 


5.2 


0.2 


0.285 


0.015 


0.0299 


0.003 


B4 


OV 


Y 


OPAL 


ro2(n) 


0.925 


0.005 


1.98 


0.07 


5.2 


0.2 


0.285 


0.015 


0.0299 


0.003 


A5 


MLT 


Y 


OPAL 


M var , A04 mix 


1.109 


0.006 


1.74 


0.06 


5.9 


0.3 


0.280 


0.013 


0.0239 


0.002 


B5 


MLT 


Y 


OPAL 


rca(n) 


0.927 


0.005 


1.84 


0.07 


5.9 


0.3 


0.280 


0.013 


0.0239 


0.002 


Ans 


MLT 


Y 


OPAL 


M var , 


1.105 


0.007 


2.40 


0.33 


8.9 


1.8 


0.259 


0.021 


0.0340 


0.003 


Bns 


MLT 


Y 


OPAL 


No Seismo 


0.934 


0.006 


2.61 


0.31 


8.9 


1.8 


0.259 


0.021 


0.0340 


0.003 



also their value for the mixing-length parameter of orCen B is 
~ 10% larger than a^. 

This agreement is not surprising, as similar observational 
constraints were considered. It strengthens the results obtained, 
since a different calibration procedure and different stellar evo- 
lution codes (with different equation of state, treatment of dif- 
fusion and treatment of sub-photospheric boundary conditions) 
were used, and it provides us with a good reference model to 
study the dependence of the calibrated stellar parameters on the 
choice of observable quantities and of the physics. 

The observation al values of the masses given by 
IPourbaix et al.l J2002I) . though precisely determined, should be 
treated as observables and therefore introduced, with their er- 
ror bars, in the definition of the x 2 and allowed to be changed 
during the calibration. Ma and Mb are considered both as pa- 
rameters and observables in a second set of fittings (A2,B2; 
A2f,B2f). The readjustment of parameters leads to a decrease 
of Mb which is 1.5cr smaller than the value determined by 



IPourbaix et alJ J2002I) . and to a decrease of age (r = 6.4 ± 
0.6 Gyr instead of 6.8 + 0.5 Gyr). The location of orCen B 
in the HR diagram has significantly improved compared with 
(Alr,Blr) (Fig.|2j and (Av A ) is also better reproduced (Fig.|3J, 
leading to an overall lower x\ compared with the equivalent 
fitting with fixed masses. 

Even including the stellar masses among the parameters, 
we are not able to improve significantly the fit of radii and large 
separations. Actually the large separation is strongly dependent 
on radius (Av oc (M/R 3 ) 1 ^ 2 ), and, given the high precision of ra- 
dius data, the procedure privileges sets of parameters providing 
the radii within lcr (for A), and 1.5<x (for B) to detriment of a 
too high large separation: Av A is always around 106.6 /uHz in- 
stead of 105.5 yuHz. In order to relax the constraints, we have 
performed a fitting including Tiff's among the observables in- 
stead of the radii (Alt,Blt). This fitting provided a small x\ 
thanks to the good match of large separation values. However, 
the radii (not included in the x 1 function) are systematically 
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Table 3. Observable quantities predicted from the sets of parameters in Table[2| 



Model 


M/M 


2 

Xri 


(Z/X) s 


2 

Xri 


R/Ro 


xl, 




7 
XRi 


L/L 




Av 




6v 




xl 


Alt 
Bit 


1.105 
0.934 




0.037 
0.039 


0.06 
0.00 


1.234 
0.874 




5764 
5280 


0.17 
0.03 


1.508 
0.534 


0.04 
0.44 


105.7 
161.1 


0.01 
0. 


4.95 
9.22 


0.66 
0.1 


1.52 


Air 
Blr 


1.105 
0.934 




0.038 
0.041 


0.02 
0.00 


1.227 
0.872 


0.20 
0.67 


5770 
5289 




1.498 
0.533 


0.13 
0.48 


106.5 
161.6 


0.52 
0.11 


5.07 
9.33 


0.53 
0.14 


2.80 


Alnd 
Bind 


1.105 
0.934 




0.040 
0.040 


0.00 
0.00 


1.227 
0.872 


0.17 
0.66 


5780 
5258 




1.508 
0.521 


0.05 
0.17 


106.6 
161.7 


0.56 
0.13 


5.34 
9.60 


0.27 
0.24 


2.27 


Ale 
Ble 


1.105 
0.934 




0.038 
0.041 


0.02 
0.01 


1.228 
0.872 


0.24 
0.66 


5772 
5286 




1.500 
0.533 


0.10 
0.44 


106.5 
161.6 


0.46 
0.44 


5.16 
9.42 


0.43 
0.17 


2.64 


Al-3 
Bl-3 


1.105 
0.934 




0.035 
0.038 


0.07 
0.01 


1.224 
0.871 


0.00 
0.29 


5767 
5315 




1.488 
0.543 


0.16 
0.51 


106.5 
161.6 


0.04 


5.77 
9.85 


0.99 


2.08 


A2 
B2 


1.104 
0.926 


0.00 
0.32 


0.039 
0.042 


0.01 
0.01 


1.226 
0.870 


0.13 
0.38 


5782 
5259 




1.509 
0.520 


0.04 
0.12 


106.6 
161.6 


0.55 
0.09 


5.35 
9.72 


0.26 
0.30 


2.20 


A2c 
B2c 


1.107 
0.919 


0.01 
0.99 


0.036 
0.040 


0.07 
0.00 


1.227 
0.870 


0.22 
0.37 


5778 
5210 




1.508 
0.501 


0.04 
0.00 


106.6 
161.0 


0.45 
0.00 


6.28 
10.54 


0.03 
0.69 


2.89 


A2f 
B2f 


1.105 
0.927 


0.00 
0.28 


0.039 
0.041 


0.01 
0.01 


1.226 
0.869 


0.12 
0.28 


5784 
5261 




1.510 
0.519 


0.03 
0.13 


106.4 
161.6 


0.32 
0.07 


5.35 
9.72 


0.27 
0.29 


1.80 


A3 
B3 


1.110 
0.927 


0.06 
0.18 


0.038 
0.042 


0.01 
0.02 


1.224 
0.869 


0.00 
0.16 


5791 
5259 




1.512 
0.518 


0.01 
0.07 


106.6 
161.4 


0.03 


5.87 
10.16 


0.94 


1.46 


A3c 
B3c 


1.113 
0.922 


0.17 
0.40 


0.038 
0.041 


0.02 
0.00 


1.223 
0.870 


0.01 
0.16 


5819 
5206 




1.540 
0.497 


0.04 
0.01 


106.8 
161.1 


0.00 


5.77 
10.22 


0.87 


1.67 


A3a G 
B3ff G 


1.114 
0.921 


0.15 
0.47 


0.038 
0.041 


0.02 
0.01 


1.224 
0.869 


0.00 
0.15 


5807 
5190 




1.529 
0.491 


0.00 
0.04 


106.8 
160.7 


0.02 


5.94 
10.35 


0.75 


1.60 


A3f 
B3f 


1.110 
0.927 


0.06 
0.16 


0.038 
0.042 


0.01 
0.01 


1.224 
0.868 


0.00 
0.12 


5790 
5260 




1.511 
0.518 


0.02 
0.07 


106.5 
161.4 


0.02 


5.85 
10.14 


0.92 


1.40 


A3nd 
B3nd 


1.108 
0.929 


0.02 
0.10 


0.040 
0.040 


0.00 
0.00 


1.224 
0.869 


0.00 
0.19 


5795 
5240 




1.516 
0.511 


0.01 
0.02 


106.5 
161.6 


0.04 


5.90 
10.16 


0.95 


1.31 


A3e 
B3e 


1.109 
0.927 


0.05 
0.16 


0.038 
0.042 


0.01 
0.01 


1.224 
0.869 


0.00 
0.17 


5792 
5256 




1.514 
0.517 


0.01 
0.06 


106.6 
161.5 


0.03 


5.89 
10.17 


0.95 


1.45 


A4 
B4 


1.112 
0.925 


0.13 
0.27 


0.036 
0.039 


0.05 
0.0 


1.224 
0.869 


0.00 
0.15 


5782 
5275 




1.503 
0.524 


0.05 
0.13 


106.7 
161.5 


0.02 


5.99 
10.56 


1.16 


1.96 


A5 
B5 


1.109 
0.927 


0.04 
0.17 


0.028 
0.030 


0.02 
0.01 


1.224 
0.869 


0.00 
0.20 


5791 
5251 




1.513 
0.516 


0.01 
0.05 


106.7 
161.6 


0.05 


5.86 
10.12 


0.95 


1.51 


Ans 
Bns 


1.105 
0.934 




0.039 
0.041 




1.224 
0.863 




5800 
5238 




1.521 
0.502 




107.0 
164.3 




3.57 
8.27 




0.06 



larger (by more than 2cr) than lKervella et all J2003ll ones. The 
large separation is very much affected by the external layers 
properties, such as either the description of the super-adiabatic 
region in the upper boundary of the convective zone, or the non- 
adiabatic processes (not taken into account either in the stellar 
models or in the oscillation code). An inspection of frequencies 
predicted for these models (Fig. [6} shows that even if the fit of 
Ava is almost perfect, the frequencies sho w a shift of 25 uHz 
with r espect to the values determined by iBouchv & CarrieJ 
ll2002l) . On the other hand, sets of parameters with a better fit 
of the radii reduce significantly the shift of frequencies. 

The small separation for the A component (Fig. 0) sug- 
gest that seismic observables would privilege younger mod- 
els, whereas one of the two values of 5v^ (n — 23) and the 
classical observables tend to a high value of the age. In fact, 
in our calibration (Ans, Bns) where only classical observable 
have been taken into account for the fit, we obtain a very good 
agreement for masses, radii, luminosity and effective tempera- 
ture (not taken as observable) for both components, and the age 
isr = 8.9 + 1.8 Gyr. 



Though the small separation averaged value of component 
A is well reproduced by our models, the slope of 5v as a func- 
tion of v differs from the theoretical prediction. One could ar- 
gue that this could b e a consequence of taking as obs ervable the 
average value of 5v. Roxbur gh & Vorontsovl ( 120031 and refer- 
ences therein), show that the Tassoul asymptotic result gives a 
poor fit both to the small separations of stellar models and to 
the observed values for the Sun, and that a better fit is obtained 
by using the ratios of small and the large separations (ro2(«))- 
This ratio depends only on the inner phase shifts which are 
determined solely by the interior structure of the star and are 
uninfluenced by the unknown structure of the outer layers. 

We have performed a similar fitting using six observational 
values of the ratio ro2 as seismic constraints for the A compo- 
nent. Unfortunately, the p-modes identified for B component 
do not allow us to define any value of ro2, we decided, there- 
fore, to take as seismic constraint the average large separation. 
The first thing to be noticed, when comparing the new set of 
parameters (A3,B3) with the one based on average large and 
small separation for both components (A2,B2), is a difference 
in the resulting age of about 1 Gyr. We also see that, by fitting 
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Fig. 3. Large (upper panels) and small (lower panels) separations for 
the A (left) and B (right) components of aCen system. Points repre- 
sent the observational values with their error bars assuming an error 
in frequencies equal to 2<r. All these curves correspond to calibrations 
performed including in the x 1 functional the average of the observa- 
tional large and small separations, but different classical observational 
constraints: masses fixed to the value given by Pourbaix et al. 1 2002) 
and r elI (solid lines); radii (dotted-lines); and masses variable as pa- 
rameter and radii (dashed-lines). All the curves were computed from 
stellar models including gravitational settling, and MLT treatment of 
convection with two mixing length adjustable parameters. For clarity 
only t = 1 theoretical large separations are shown in the upper panels. 

ro2, we get a good fit of 6vo2 for both stars (Fig.@J. The large 
separation is slightly higher than the observational one, as we 
obtained in (A2,B2) and (Alr,Blr). 

In Fig. [6] we can see also the effect on the p-mode frequen- 
cies, the model (A3,B3) providing better fit to the observational 
ones, than the model (A2,B2). 

Notice that the parameter set (A3,B3) was determined with- 
out taking into account the small separation for the B compo- 
nent. Theoretical (6v) is ~ 10.15 uHz, that is close to the new 
observational value (10.1;uHz) bv )Kjeldsen & Bedding (2004). 
We wonder whether the high age obtained by takin g (5vb) as 
given by the average of two points JCarrier & BourbaiJ2003l) is 
only a consequence of the low value imposed to (dv^). In fact, 
a different calibration performed using the small and large sep- 
arations as observables, but taking only 6v^{n = 21) instead of 
the average of two available values, provides an age of 6.0 Gyr, 
in good agreement with ro2 (A3,B3) fittings, and quite younger 
than (A2,B2) models. 

Finally, we have also used ro2 but without varying the 
masses. Again, the age is of the order of 6.0 Gyr, but now, 
the mixing-length parameters needed for both stars are quite 
different (by 17%). 

The fits reported in Table[3]could indicate an inconsistency 
between the seismic and classical observables. Actually for 
component B the resulting radii are larger than the observed 
ones by more than ~ 1 cr or ~ 2<x and the masses are smaller by 
more than ~ 1 cr or ~ 2cr with respect to the values observation- 
ally derived. These results could be interpreted as a systematic 
error in the radii and/or in the mass determination. A calibra- 



Fig. 4. As Fig.|3|but for models computed using different approach for 
convection: MLT with two mixing length parameters (solid line); FST 
(dashed lines); MLT with a A = a B (dotted lines); and MLT with a A = 
= a Q (dash-dotted lines) Large and small separations are derived 
from stellar models, but the fitting is performed using the ratios r 02 for 
component A and Av for component B. 



tion performed with a free Mb parameter/observable, leads to 
M B = 0.91 1 M and R B = 0.863 R . The other parameters do 
not significantly change. This mass value would imply a sys- 
tematic error of ~ 4<x that, given the available high precision 
data, does not seem reliable. On the other hand, if R% is let 
free, the fit of the other observables provides R B = 0.871 R 
instead of O.863R , and a M B value within 1 cr from the ob- 
served one. It must be notice, however, that these fits have in- 
cluded Avb in the^- 2 functional. To check if the large radius is 
only a consequence of Avb, we calibrate the system adopting 
as seismic constraints ro2(«) for component A, and 6v(n = 21) 
for component B. In this case we are able to fit the radii and 
masses of both components within lcr, but (Avb) is more than 
+3-cr far away from the observational value, and the predicted 
frequencies are ~ 40 /jHz larger than observed ones. A pri- 
ori, we cannot rule out a larger uncertainty in the observed ra- 
dius. Actually in the observational radius determination there 
are implicit a definition of stellar radius (that is not necessarily 
the same than that used in stellar modelling) and an assump- 
tion about the limb darkening law and the atmosphere models. 
How much the stellar radius is affected by these assumptions? 
Kervella (2005, private communication) claims that these ef- 
fects are much smaller that other errors intrinsic to the mea- 
surement method and already taken into account. So, we should 
wait for more precise seismic data for o/CenB to understand this 
apparent inconsistency. 

5.2. Effect of the treatment of convection 

It is well known that one of the weak points of stellar evolution 
models is the treatment of convection. The "standard model" of 
convection adopted in stellar evolution is the mixing length the- 
ory (MLT) where turbulence is described by a relatively sim- 
ple model that contains essentially one adjustable parameter: 
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the mixing length A = aH p (H p being the local pressure scale 
height and a an unconstrained parameter). The value of this pa- 
rameter determines the radius of the star and the behavior of the 
super-adiabatic region in the outer boundary of the stellar con- 
vective zone. The calibration of a stellar code using the solar 
radius and solar luminosity provides the value of the mixing- 
length parameter (a Q ) which in turn is usually used to model 
other stars with the same physics. A question arises: can stars of 
different masses, initial chemical composition and evolutionary 
status be modeled with a unique al If the answer is negative, 
how reliable is the shape of the isochrones determined with a 
unique al 

o-Cen offers a unique opportunity of testing our as- 
sumption about a and, therefore, our simplified way of 
overcoming the complex problem of convection in stars. 
The question of the universality of the mixing-length pa- 
rameter has been approached many t imes in the past 
lNoelset al]ll99l| lEdmonds et al.1 Il992t iLvdon et all Il993t 
Neuforge 1993: iMoreletalJ 2000: iFernandes & Neuforg d 
1995; Guenther & Demaraua l2000n . Some attempts of cali- 
brating aCen A&B in luminosity and radii using MLT theory 
suggest different values of a for each component and differ- 
ent from the a Q , while others favor similar values for both 
components, or conclude that the uncertainties in masses and 
radii as well as the chemical composition should be reduced 
significantly before being able to draw fir m a conclusion on 
whether the MLT parameter is unique or not JLvdon et al.ll 993 ; 
lAndersenlll99lt iGuenther & DemaraudbooaT Nowadays, the 
high precision with which we know masses and radii for 
aCen A&B allows to analyze again this problem. The frequen- 
cies are very sensitive to R and therefore to a, and the differ- 
ence between the observed frequencies (v bs) and the theoret- 
ical ones (vtheo) 1S highly affected by how the super-adiabatic 
zone is described (see e.g. lSchlattl et al.ll997t) . 

In order to analyze these questions we have made several 
fits of the aCen observable quantities by using MLT with i) 
a a + ffB as free parameters; ii) a\ = qtr free parameter ; and 
Hi) oa — «b = «o fixed; and iv) since Mo rel et alJ (12000) pre- 
sented also controversial results when comparing calibrations 
for QfCen using MLT or FST, we shall analyze as well the effect 
of using the FST treatment of convection. 

5.2.1. MLT versus FST 

Mimicking the spectral distribution of eddies by one "average" 
eddy, such as MLT theory does, has crit i cal co nsequences on 
the computation of MLT fluxes. Canuto (1996) shows that in 
the limit of highly efficient convection MLT underestimates the 
convective flux , and on the contrary, in the low efficiency limit, 
MLT overestimates the convective flux. FST models attempt to 
overcome the one-eddy approximation by using a turbulence 
model to compute the full spectrum of a turbulent convective 
flow. As consequence, the convective flux is ~10 times larger 
than the MLT one in the limit of high efficiency, and ~0.3 in 
the limit of low efficiency. This behavior yields, in the super- 
adiabatic region at the top of a convective zone, steeper tem- 
perature gradients for FST than for solar-tuned MLT. The mix- 



ing length used in FST treatment is in general defined as the 
distance from a given point to the boundary of the convective 
zone. However, the FST fluxes have also been used combined 
with o ther definitions of mixing length. For instance, Bernkopf 
use t he FST fluxes with a mixing length A = a*H p 
with a* < 1 . iMorel et"afl d2000l) use also this kind of descrip- 
tion of FST convection with CM91 for the fluxes. Here we will 
use this prescription of FST treatment of convection, with the 
fluxes bv lCanutoetal](ll996l) instead of CM(91,92). The main 
difference is in the temperature gradient, steeper in CM than 
in CGM9 6, but both, in any case, much steeper than the MLT 
one. lMorel et al-l lEoOO) found significant differences in the fit- 
ted parameters (in particular a 1.5 Gyr difference in the age) 
when investigating the effect of using a different theory of con- 
vection in the modelling of a Centauri. This is rather surprising 
since the treatment of convection affects only the very external 
layers of the star, and such an important effect on the age of the 
star is not expected. 

The convection parameters that result from our fits are, as 
in the case of the calibration with MLT, close to the parameter 
needed to calibrate the Sun (0.753). Furthermore, in compar- 
ing the models (A3,B3) with (A3f,B3f) ones, we do not ob- 
serve any difference in the parameters of the models fitting the 
system orCen using MLT or FST. In principle this is a logic 
behavi or since we considered as s eismological observable the 
ro2(n) (Roxburgh & Vorontsov 2003) ratio, and this quantity is 
independent of surface properties of the model. However, also 
the sets of parameters (A2,B2) and (A2f,B2f) determined by 
using large and small frequency separations as constraints, are 
in fact the same, independently of the convection description 
used. 

We would expect, like in the Sun, a decrease of the differ- 
ence between theoretical and observational frequencies in the 
high frequency domain when FST is used. However, for our bi- 
nary system, the improvement provided by FST treatment with 
respect to MLT is quite small. 

For aCen A the improvement is of the order of 4 pMz at 
3000 yuHz, while for the B component, this effect does not ap- 
pear clearly (Fig|6j. In fact, the difference between MLT and 
FST frequencies is due to a different radius and mass. If we 
adopt for Mb and Rb the values determined in the calibration 
B3 as the real ones, and compute a new FST calibration, we 
find that the difference between the MLT and FST frequencies 
in the observational domain (3000-4600/iHz) goes from to 
3 yuHz, while the difference in frequencies between B3 and B3f 
introduced by the different mass and radius is of the order of 
8 //Hz. 

5.2.2. One or two different values for the mixing 
length parameter? 

Several times in the literature it has been addressed the ques- 
tion if convection in the two components of a Centauri 
should be described with distinct mixing-length parameters 
and, if, given the observational uncertainties, the inferred dif- 
fer ence between mixing leng th parameters is significant (see 
e.g lEggenber ger et alJ (120041) . where an exhaustive review on 
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previous calibratio ns is also presented). It was, in fact, already 
suggested in Guenther & Demarque (2000) that a reduction in 
the observational uncertainties and the inclusion of seismic 
constraints in the modelling would allow a more robust infer- 
ence on the mixing-length parameters of ffCen A and B. This is 
now the case, thanks to the detection of solar-like oscillations 
and to the precise determination of radii. These are compatible 
with effective temperatures spectroscopically determined and 
significantly reduce the error box in the HR diagram. As a gen- 
eral result of the calibrations presented in the previous sections 
we find that the mixing-length parameter of calibrated model 
A (ota) is approximatively 5-10% smaller than a^. 

Numerical simulations of convection has recently permitted 
to carry out a calibration of mixing-length parameter through 
the HR diagram. A function a(T e fi , log g) is determined in order 
to reproduce the step of the specific entropy provided by the 
atmosphe re hydrodyna mic models. Both calibrations, the one 
by Ludwig et al. ( 1999) (based on 2-D simulations) and the one 
by Trampedach (2004) (based on 3-D simulations) show slight 
variations of a with the position in the HR diagram, and suggest 
that the mixing parameter should be represented as a function 
of log r e ff, log g, and chemical composition. 

Actually, comparing our results with these theoretical pre- 
dictions, we see that the difference between the value deter- 
mined for an and a a is in good agreement with the predictions 
by numerical simulations of convection. Moreover, our a values 
for the two components bracket that obtained by calibrating the 
Sun (a G = 1.91) with same physics (see for instance the model 
A3,B3). This is what one expects from their HR location and 
from the calibrations by Ludwi g et alJ II1999B and Trampedach 
We note, however, that in the fits obtained by using 
(Av) and (5v) we find a% > but also larger than the 
solar one. The difference with respect to the solar value is even 
larger when the masses of the com ponents are assumed to b e 
fixed. The same effect is present in Eggenbergeret Hli004): 
they find a% > oa, but their values are far from their solar one. 

In order to determine whether the addition of an extra free 
parameter leads to a significant improvement of the fit, we per- 
formed a calibration assuming a single mixing-length parame- 
ter for both components, and including the ro2(«) ratios in the 
quality function (A3c,B3c). The value of a for this calibration 
is a = 1 .95 (quite close to the solar one), but the masses now 
are slightly larger for the component A, and slightly smaller 
for the component B. The same happens if we decide to fix 
a a = as =0(0. The masses are anyway within 2<x of the ob- 
served value, and the age is the same r = 5.7 - 5.8 Gyr. 

A different result is obtained if the small and large differ- 
ences are considered in the quality function (A2c,B2c), and the 
mixing-length parameter is assumed to be the same for both 
stars. In this case, the value reached (a = 1.78) is not so close 
to the solar one, Mb must decrease to O.919M , and the age of 
the system decreases with respect to the value obtained allow- 
ing two differen t mixing-length parameters . This result recalls 
that obtained bv lThevenin et alJ d2002h . who using a unique a 
had to decrease the M B to 0.907M o . 

Since adding a free parameter would naturally lead to a bet- 
ter fit (a lower y 1 as defined in Eq.Q), for a quantitative com- 
parison between the quality of the fit obtained with a different 
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Fig. 5. As Fig. but with different curves corresponding to differ- 
ent physics included in the stellar comput ation: convective overshoot- 
ing (dashed lines); solar mixture from Asplund et al. 1 2004) instead of 
iGrevesse & NoelJ <1993t) (dash -dotted lines); no gravitational settling 
(solid lines); CEFF equation of state instead of OPAL01 (dotted lines). 



number of model parameters, it is more meaningful to compare 
the so-called "reduced"^ (Eq.©. As can be seen in Table|3] 
the value of x\ is lower if two distinct mixing length param- 
eters are used in the modelling (both comparing A3,B3 with 
A3c,B3c and A2,B2 with A2c,B2c): this suggests that, with the 
adopted observational constraints, the addition of another free 
mixing-length parameter is justified. In fact, comparing the fits 
obtained with one or two mixing-length parameters, the F-test 
gives a confidence of 85-90% that the inclusion of two differ- 
ent parameters for convection is significant. We should how- 
ever recall that such a statistical test, and generally a y 2 statis- 
tics, assumes that the observational errors are distributed about 
the mean following a Gaussian distribution. This is not neces- 
sarily true as systematic shifts in the observed quantities and 
inaccuracies in the models cannot be excluded. As could also 
be expected, the uncertainties adopted with the observational 
constraints are crucial. For instance, we find that if the obser- 
vational error in the mass of the component B is doubled, the 
addition of a second free parameter for convection is no longer 
justified. 



5.3. Disentangling the physics 

One of the principal motivations for pursuing the study of pul- 
sations in other stars is to test the assumptions concerning the 
physics underlying the stellar structure theory. To practically 
investigate this idea, we compute four models of aCen that 
use the same observational constraints and model parameters 
as the reference calibrations (A3,B3) and (Alr,Blr), but incor- 
porate changes in the phy sical assumptions. Th ese are : CEFF 
EoS (Ale,Ble; A3e,B3e), Asplund et all ( 12004 chemical com- 
position (A5,B5). To test the effects of gravitational settling 
we have calibrated the models (Alnd,Blnd) and (A3nd,B3nd) 
without microscopic diffusion. 
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Fig. 7. Difference of frequencies between a reference calibra- 
tion (A3,B3), and both, that with a different EoS (A3e,B3e) 
(solid lines), and without microscopic diffussion (A3nd,B3nd) 
(dahsed lines). Upper panel corresponds to component A, and 
lower panel to component B. 



Fig. 6. Difference between theoretical and observed frequen- 
cies for the sets of parameters whose separations have been 
plotted in Fig. [3] (upper panel); in Fig. |4] (middle panel); and 
Fig ID (lower panel). Left side corresponding to aCen A, and 
right side to aCen B. The labels correspond to the identifica- 
tion of the model used in Table|2] 

The central question is whether the residual errors are ob- 
servationally significant. If the residuals are all small com- 
pared to the observational errors, it is not possible to distin- 
guish parameter changes from changes in physical assump- 
tions. In general, we find that stars obeying different physics 
succeed remarkably well i n masquerad ing as stars that merely 
have different parameters (Br own et all 1 9941) . 

5.3.1 . Equation of state 

Concerning the equation of state, the calibration assuming 
CEFF (A3e,B3e) leads to a result in agreement with the fit 
(A3,B3) computed adopting OPAL01, both in terms of ^ (ei- 
ther the total one or that corresponding at each observable 
quantity), and of the fitted parameters. This is expected for 
stars with an int ernal structure similar to the sun: as shown 
bv lMighoM2004 the differences between the sound speed and 
the first adiabatic exponent (Fi), due solely to the use of a 
different EOS (CEFF and OPAL01), are smaller than 1% ex- 
cept in outer regions with radii larger than 0.95 R«, that is of 
the same order of those predicted i n solar models (see e.g. 
iBasu & Christensen-Dalsgaardifl99 7). The differences appear 



to be larger in the lower-mass model B than in model A, and are 
mainly located in the hydrogen and helium ionization regions. 
In that study, the differences in T\ came from the application of 
CEFF or OPAL01 to a given stellar structure (logp vs. log T) 
and chemical composition. In the calibration procedure, such 
small differences in the internal structure of a model propagate 
in a variation of the observables of each model that could be 
easily compensated by a re-adjustment of the free parameters. 

The models calibrated by using different equations of state 
are very similar. There is, nevertheless, a difference in the fre- 
quencies. In Fig.[7](solid lines) we plot the difference between 
(A3,B3) and (A3e,B3e) frequencies. The upper panel corre- 
sponds to the A component, and the lower panel to the B 
one. This difference in frequencies shows a behavior depend- 
ing on v with an oscillatory signature. This is expected for in- 
stance, when the models have either different locations of the 
convective region boundaries, or a different behaviour of T\. 
The bottom of the convective zone in o-Cen A is located at 
*zc - Rcz/R* = 0.707 for A3 calibration, and x zc = 0.711 for 
A3e model. Therefore, the oscillatory signature is probably re- 
lated to the changes in the second helium ionization zone. The 
amplitude of this oscillation is linked to the difference in depth 
of the F[ bump, and the period contains information about the 
location of the difference in stellar structure. For star A, the 
amplitude of the oscillation between the order n — 10 and 15 
is around 0.7 /jHz, while for the B component that is almost 
2 //Hz. 
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Fig. 8. no()i) ratios for A component. Points represent the observa- 
tional values with their error bars assuming an error in frequencies 
equal to 2a. Short-dashed line corresponds to the model calibrated 
using large and small separation in the x 1 functional, as well as the 
effective temperatures instead of the radii. The other three curves 
correspond to modes calibrated using r 02 (n) in x 1 '- model including 
overshooting (solid line) and presenting a convective core; the same 
calibration without overshooting (dot-dashed line), and calibration 
without overshooting and without microscopic diffusion (long-dashed 
line). 

5.3.2. Microscopic diffusion 

The calibrations without diffusion (Alnd,Blnd) and 
(A3nd,B3nd) provide fits of similar quality with respect 
to the corresponding calibrations including gravitational 
settling, respectively (Alr,Blr) and (A3,B3). As reported in 
Table [2] the parameters resulting from the calibrations differ, 
as expected, in the initial chemical composition. The smaller 
mass fraction of He required in these calibrations implies 
also a slight increase in the age of the system: 6.3 + 0.4 Gyr 
for (A3nd,B3nd) versus 5.8 + 0.2 Gyr for (A3.B3) (and 
7.1 + 0.6 Gyr for (Alnd,Blnd) versus 6.8 + 0.5 Gyr for 
(Alr,Blr)). The (A 3nd,B3nd) c alibration is consistent with 
the one obtained bv iThoul et alJ J2003h . Since they use only 
cCen A frequencies, their age estimation is not affected by the 
low 6v(B) value. 

The mixing-length parameters are also re-adjusted to fit the 
stellar radius, and also in this case is almost 8% larger than 
qta, and both bracket the mixing length parameter obtained for 
a Sun calibrated without microscopic diffusion (a e = 1.84). 

In Fig. (dashed-line) we show how diffusion affects 
the frequencies. We plot the residuals defined as v„ (A3) - 
v„o(A3nd) (upper panel) and v„ (B3)-v„o(B3nd) (lower panel). 
For component A the residuals as function of the order n show 
an oscillatory behavior, reflecting the changes in envelope He 
abundance (F s ). Actually the calibrated model including dif- 
fusion has a superficial He abundance F S (A3) = 0.245, while 
model A3nd, without gravitational settling has Y s = 0.270. This 
difference implies as well a different opacity and, therefore, 
a different location of the boundary of the convective zone 
(X cz (A3nd) = 0.725). To avoid changing the figure scale, the 



curve of residuals in the upper panel has been shifted down by 
2.5^Hz; we see that the amplitude of the oscillatory signal is 
around 0.5 /iHz between n = 10 and n = 20. 

For component B the curve of residuals shows a com- 
pletely different behavior. Since orCen B is less massive than 
its companion, the mass contained in its convective zone is 
much larger and, therefore, the effect of microscopic diffusion 
is much smaller. The value of Y s for B3 is less than 4% smaller 
than that for B3nd, while for the component A, the difference 
in F s between both calibrations (A3 and A3nd) is around 10%. 

5.3.3. Solar mixture 

Asol und et all (12004 have recently proposed a substantial re- 
vision of the abundance of C N and O in the solar photosphere. 
Since the value of the surface metallicity of o-Cen is observa- 
tionally determined relative to the solar heavy element abun- 
dance, the solar (Z/X) s = 0.0177 resulting from the new so- 
lar calibration would lower (Z/X) s in aCen by ~ 30%. We 
thus calibrate our models assuming as observational constraint 
(Z/X) s = 0.029 (instead of 0.039) and using in the computa- 
tions OPAL opaci ty tables c alcu lated with the solar mixture 
proposed by Asp lund et alJ J2004I) . The results of this cali- 
brations, models (A5,B5), present a different initial chemical 
composition, but no other significant deviation from the global 
parameters of models (A3,B3) is noticed. As it happens for 
the Su n, the aCenA model calibrated with iGrevesse & Noelsl 
(1993) has a deeper convective region compared with the model 
A5, but the uncertainty in the observational data is quite larger 
than in the Sun, and this difference can be masked by other 
choice of parameters. 

5.3.4. Overshooting 

An additional consequence of our simple way of describing 
convection in stellar models, is the need to param eterize con- 
vective overshooting. In gener al (see e.g. lSchaller etalJI 19921) 
convective core overshooting is necessary to fit isochrones of 
open clusters for masses larger than a given critical mass. The 
problem is to determine this critical mass and to describe the 
transition between no-overshooting mass domain and over- 
shooting mass domain. On the one hand, the mass and ef- 
fective temperature of aCen A are quite close to solar val- 
ues, and one could think that no overshooting should be in- 
troduced. Nevertheless the chemical composition of o-Cen A 
is different from the solar one and the evolution of a convec- 
tive core does not necessarily show the same behaviour. The 
mass of aCen places it in the boundary region between mod- 
els with and without a convective core. We have made sev- 
eral calibrations varying the thickness of the overshooting layer 
ov — p * (min(r cv , H p (r cv ))) (where r cv is the radius of the 
convective core) with p from 0.0, 0.1, 0.15 and 0.2. In fact, 
for values of p < 0.15 no convective core remains after the 
PMS, therefore the parameters resulting from the fitting are not 
changed. 

In Tables |2] and [5] we report the set of parameters and 
observables corresponding to models (A4,B4) calibrated with 
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P = 0.2. As should be expected, including overshooting re- 
duces the age of the system. Fig. [5] shows the large and small 
separations. We find a possible direct indicator of a convective 
core in the behavior of the small separation of component A 
at high frequencies (v > 2500 yuHz). That is even more ev- 
ident in the signature left in ro2- However, given the uncer- 
tainty (1.3/jHz) affecting the p-modes involved in the highest 
frequency ro2(«) (or Sv{n)), it is not possible, based only on 
ro2(«) or c)vo2, to rule out a convective core in aCen A. The 
ratio no(n) = (v„o - 2v„i + v„+i,o)/(2Avo(n + 1)), however, is 
much more eloquent, as shown in Fig.|S] and allows us to reject 
model A4: current observational constraints seem not to be in 
favour of a convective core in aCen A. 

6. Summary and conclusions 

We propose a calibration of the binary system aCen by means 
of Levenberg-Marquardt minimization algorithm applying it 
simultaneously to classical (photometric, spectroscopic and 
astrometric) and seismic observables. The main features of 
this sort of algorithm make it an ideal tool for the aim of 
this work: to practically analyze the effectiveness of oscilla- 
tion frequencies in constraining stellar model parameters and 
stellar evoluti on physics , by using the p- modes identified by 
iBouchv & Carrierl d2002t) (aCen A) and bv lCarrier & Bou rban 
fcOO^ faCen B). Actually thanks to its low computational cost, 
this algorithm allows to search for the best model in the full 
7-dimensional parameter space describing the binary system, 
varying both the set of stellar observables, and the physics in- 
cluded in stellar modelling. 

As starting po i nt we assumed the same observables than 
Eggenber ger et alJ J2004. In spite of the different EoS (MHD), 
diffusion treatment ( Ri chard et all 19961 for five elements sep- 
arately), cross section of nuclear reactions (NACRE) and atmo- 
spheric boundary conditions, we derive a set of stellar parame- 
ters (Alr,Blr) in complete agreement with theirs. 

By comparison between calibrations with different observ- 
ables we make clear that care has to be taken when using Avo 
to constrain fundamental stellar parameters. Given the strong 
dependence of Avo on surface layers, and our poor understand- 
ing of the physics describing outer stellar regions, seeking a 
perfect agreement between observed and predicted Avo may 
bias our results toward, for instance, inaccurate radii. This is 
clear when comparing the models Alt and Air; in the former a 
perfect agreement with the observed Av is reached but the radii 
(and the frequencies, see Fig.|6j deviate significantly from their 
observational values. 

An additional source of systematic error in the calibrations 
concerns the value of the small frequency separation of com- 
ponent B. Since the observational data are still rather poor one 
of the two measured values, 5v^{n - 23) leads to a higher age 
than suggested by 6va- 

The value of Sv% predicted in our calibration (A3,B3) 
(where Sv% is not included in the x 2 ) is m agreement with 
6vn(n - 21) and with th e very recent value published by 
iKieldsen & Bedding| (Eo04). This also suggests that sufficiently 
precise seismic data of one star are sufficient to determine fun- 
damental parameters of the system. In fact additional calibra- 



tions, not shown in Table |2 where no seismic constraints of 
component B are included, lead to fitted parameters compati- 
ble with e.g. (A3,B3) even though a large discrepancy between 
predicted and observed frequencies of aCen B is observed. 

We therefore propose to use ro2 as a reliable seismic con- 
straint to determine fundamental parameters of a star. The large 
frequency separation could, nonetheless, prov ide a first esti - 
mate of the mean density and, as shown e.g. in Gough ( 1990), 
an useful information on localized features in stellar interior 
once more accurate determinations of solar-like oscillations 
will be available. 



Section 15721 was devoted to the study of the effects on the 
calibration of our uncertainties concerning stellar convection. 
If the calibrations, e.g. (A3,B3), are performed considering a 
free parameter describi ng convection in ea ch component (ota, 
as) we find, as lE ggenberger et all (|2004), a s 5-11% higher 
than »a- Differently from Eggenberger et al. (2004) we find 
that a a < a G < a%: this is of primary relevance when mak- 
ing statements concerning the difference between a a and cb, 
otherwise we would also have to justify an even bigger (not ex- 
pected) difference between a\ and a Q . We notice also that our 
inferred values of a follow the sa me trend pred icted by MLT 
parameter calibration based on 2D jLudwigetaill9 99) and 3D 
(Trampedach 2004) hydrodynamic atmosphere models. 

In order to draw more firmed conclusions on the signifi- 
cance of considering a- a and qtb as free parameters we repeated 
our calibrations (A3,B3) and (A2,B2) assuming a single pa- 
rameter for both components ((A3c,B3c) and (A2c,B2c)). The 
fit necessarily improves when an additional free parameter is 
introduced in the calibration, nevertheless we find, with the ob- 
servational constraints we adopted, the difference between a\ 
and »b significant. We note also that the value of a obtained 
in (A3c,B3c) and (A2c,B2c) calibrations is quite close to the 
corresponding a G . 

We also find that, contrary to what was obtained by 
iMorel et"al] d2000l> . the use of a different theory of convection 
(MLT or FST) in our models does not change the set of pa- 
rameters derived from the fitting. The only effect of convection 
model is a slight improvement in the fit of high frequencies 
when using FST (in particular for A component). 

The available seismic data are not in favour of a convective 
core in a Cen A, moreover, the overshooting parameter needed 
for a convective core to persist after the PMS (fJ = 0.2) appears 
to be too large for a mod el of the mass and che mical composi- 
tion of a Cen A (see e.g. lDemarque et al.l2004l) . 

Finally, we find that stars obeying different physics provide 
similar fits to those obtained with stars that merely have differ- 
ent parameters. As a consequence, we are not able from present 
data to discriminate neither between different EoS, neither be- 
tween diffusion or not diffusion models. Fig. shows, how- 
ever, that expected precision from space missions would allow 
to apply inversion techniques to frequency data. Adding seis- 
mic observables has significantly improved the determination 
of the system fundamental parameters, but more and more pre- 
cise observations are needed to be able to extract information 
about the internal structure of the stars. 
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